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Abstract 

We consider a one-dimensional system consisting of a granular tracer particle of mass 
M in a bath of thermalized particles each of mass m. When the mass ratio, M/m, 
is equal to the coefficient of restitution, a, the system maps to a a one-dimensional 
elastic gas. In this case, Boltzmann equation can be solved exactly. We also obtain 
expressions for the velocity autocorrelation function and the diffusion coefficient. 
Numerical simulations of the Boltzmann equation are performed for M/m ^ a 
where no analytical solution is available. It appears that the dynamical features 
remain qualitatively similar to those found in the exactly solvable case. 
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1 Introduction 



Granular gases consist of particles that undergo dissipative collisions, making 
them a paradigm of non-equilibrium statistical mechanics 0, 0| • The absence 
of microscopic reversibility (in each collision, kinetic energy is lost) and the 
contraction of the phase space volume lead to specific properties which differ 
from elastically colliding thermal systems. These include non-gaussian statis- 
tics, modified hydrodynamics, and absence of equipartition [lfM HI 111 • 

Compared to equilibrium, our knowledge of nonequilibrium statistical mechan- 
ics is still incomplete, and exact results are rare. In this paper, we consider a 
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one-dimensional model of a tracer particle that undergoes dissipative collisions 
with particles of a thermalized bath. The kinetic description is provided by 
the Boltzmann equation. We show that when the mass ratio of the tracer to 
a bath particle, M/m, is equal to the coefficient of restitution a, there exists 
a mapping with the Boltzmann equation of identical elastic hard rods, the 
solution of which was derived almost thirty years ago 0]. We propose here 
a simpler way of solving the Boltzmann equation and we discuss the results 
in the context of our model, i.e. a granular particle in a thermalized bath. 
In addition, we obtain the velocity autocorrelation function as well as the 
diffusion coefficient. Numerical simulations of the Boltzmann equation show 
that the qualitative kinetic features are not strongly dependent on the value 
of the coefficient of restitution, implying that the exact solution provides the 
characteristics of the system also when M/m ^ a, where an analytic solution 
is not available. 



2 The model 



We consider a one-dimensional gas of identical hard rods of mass m into 
which one inserts a granular hard rod of mass M0|. Particles have hard core 
interactions and collisions between bath particles are elastic. But, collisions 
between the tracer particle and the bath particles are assumed inelastic and 
characterized by a coefficient of restitution, a < 1. In addition, we assume 
that a stationary, Gaussian velocity distribution is imposed on the bath parti- 
cles, compensating for the kinetic energy lost during collisions with the tracer 
particle. 

For each collision, momentum is conserved so that for a collision between a 
bath particle and the tracer particle, one has 



where the upper-case velocity corresponds to the tracer particle and the lower- 
case velocity to the bath particle, and the asterisks denote post-collisional 
quantities. At the moment of impact, the relative velocity changes sign and 
shrinks (for a < 1) according to the collision rule 



By combining Eqs.(l) and (2), the velocities of colliding particles after collision 



MV* + mv* = MV + mv 




V* -v* = -a(V-v) 



(2) 



2 



are given by 

v * = (fJ.-a)V+ (1 + a)v ,gx 
1 + n 

^ = (l+q)F + ( / /- 1 -q)^ 
1 + /X -1 

where /i = M/m. 

Note that the tracer mass can be also taken as e qual to the fluid particle mass 
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by introducing an effective restitution coefficient 
The precollisional velocities V** and v**, corresponding to the inverse collision 

(V**,v**) -> (V», 

can be expressed in terms of V and v by replacing in Eqs.(3) and (4) the 
coefficient a by a' 1 . 

3 Structure of the Boltzmann equation at a = M/m 



The kinetic description of this model is provided by the Boltzmann equation 
according to which the probability distribution of the tracer particle evolves 

as 

/+oo . , 

dc\V - c\ \a~ 2 f(R, V*\ t)0(c**) - f(R, V, t)0(c) (5) 

where R is the position of the tracer particle, <fi(c) denotes the equilibrium 
Maxwell distribution of the bath particles, and p their number density. Note 
that Eq. (5), although linear in the case of a tracer particle, has not yet been 
solved exacty. Exact solutions have, however, been found in the framework of 
the inelastic Maxwell model where the collision frequency, which occurs in the 
Boltzmann collisional term, is assumed independent of the relative velocity 

SE3- 

The first objective of this paper is to show that an exact solution of Eq. (5) 
can be obtained when the coefficient of restitution equals the mass ratio of the 
tracer to a bath particle, a = fi = M/m. In this case, the Boltzmann equation 
for the distribution of the granular tracer particle turns out to be mathemat- 
ically identical to the Boltzmann equation describing a tagged elastic particle 
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in a thermalized bath of mechanically identical particles. This fact, shown be- 
low, opens the way to an exact solution. However, the physical situations are 
totally different: evolution of a dissipative system towards a stationary state in 
one case and evolution of a conservative system towards thermal equilibrium 
in the other. 

The differences between the two systems on the level of the microscopic dy- 
namics are illustrated by the collision rules. For the elastic system with iden- 
tical masses, a binary collision leads to the exchange of particle velocities, 
V* = c and c* = V. Moreover, for the inverse collision one finds V** = c and 
c** = V, which is related to the microscopic reversibility. For the granular 
tracer particle of mass M = am in a thermalized bath of particles of mass m, 
the collision rules are c* = aV + (1 — a)c and V* = c, whereas for the precol- 
lisional velocities corresponding to the inverse collision one finds c** = V and 
V** = c+ (1 — a)V. Despite the differences, in the gain term of the Boltzmann 
collision operator, in both cases, the precollisional velocity of the bath particle 
equals V, and is thus independent of c. And this common feature results in 
the simplification that allows the equation to be solved, as can be seen below. 

For the sake of simplicity, we first consider the Boltzmann equation for spa- 
tially homogeneous systems. The general solution for inhomogeneous initial 
conditions can be derived in a very similar way. The corresponding, rather 
tedious, calculations are given in Appendix A. 

The Boltzmann equation for homogeneous states takes the form 



(V,t) = p dc\V-c\ 




1 + a'^V + - a 
1 + 



'(p-a- 1 )V+ (1 + a~ 1 ) 
- f(V,t)<f>(c) 



(6) 



Putting \i = a in the gain term of Eq.(6) yields 



' d_ 
dt 



+ u{V) I f(V, t) = pcj)(V) I dc\V- c\a~ 2 f I ^ ^ +-,t] (7) 




where v(V), the collision rate of the tracer, is given by 



+oo 



u(V) = p / dc\V - c\(f)(c) 



By using a change of the integration variable, u = c/a + (a — l)V/a, in the 
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right-hand side of Eq.(7), one finally obtains 



(J^ + v{v)jf(y,t) = p<f>(y) J + ™du\v-u\f(u,t), (9) 

which coincides exactly with the Boltzmann equation for elastically colliding 
identical particles of mass m, solved by Resibois 0]. 

Before determining the time evolution of f(V, t), we first focus on the station- 
ary state f st {V) which satisfies the equation obtained by dropping the time 
derivative in Eq.(9): 

/+oo r+oo 
dc\V - c|0(c) = (f)(V) / dc\V- c\f st (c), (10) 
-oo J — oo 

The solution of (10) is given by 



where T is temperature of the bath. Therefore, the stationary velocity dis- 
tribution is Maxwellian, but with a temperature aT smaller than the bath 
temperature. This result is a specific case of a previously obtained general 
formula valid for any mass ratio in any dimension [lj|. 

From Eq. (8) we get an explicit expression for the collision frequency of the 
granular particle 



v{V) = p 



2aT f MV 2 \ T/ , 
exp — — — + V erf 




:i2i 



7iM 1 \ 2aT 
where erf(x) is the error function. 

In the following, we use dimensionless variables: the units of length and time 
are 1/p and y/m/VTp, respectively. Also, for the sake of simplicity, we use 
lower-case variables since no confusion with the variables associated with the 
bath particles is possible. 



4 Exact Homogeneous solution 



We first show that the integro-differential equation (9) can be replaced by a 
differential equation. To this end we introduce an auxiliary function G(u, t) 

/+oo 
dc\u-c\f(c,t) (13) 
-oo 
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satisfying the relation 

|^(M) = 2/(M). (14) 

Therefore, the Boltzmann equation, Eq.(9), when re-expressed in terms of 
G(u, t), yields 

d 3 G , , , ,d 2 G , , d 2 v 



mdu^ t) + u ^ M -^ u)GM (15) 

Let us introduce the time Laplace transform 

r+oo 

G(u,z) = / dtG(u,t)e- zt , (16) 

JO 

which when applied to Eq.(15) yields 



, 1 d 2 G , , d 2 v{u) ~, , d 2 G . 
[Z + U( - U)] &^^ Z) ~ Z) = 0) (17) 



By using the identity 



-H|^K-)-^HG(n,,) = |- 



we can integrate Eq.(17) over the velocity from — oo to u finding 

[z + z) - ^(u)G(u, z) = ^-(u, 0) + B(z) (19) 

where B(z) is a function that can be determined by taking in (19) the limit 
u — * oo. According to the definitions (13), (16) in the region of large velocities 
G(u, z) is asymptotically given by 

~ \u\ 

Giu, z) — — - sgn(u) < viz) > (20) 

z 

where < v(z) > is the Laplace transform of the mean velocity 

/+oo 
duuf(u,z). (21) 
-oo 

Using Eq.(20) one can readily check that Eq.(19) considered in the limit |it| — > 
oo reduces to the equality B(z) =< v(z) >. 



Introducing now the function 



Giu, z) . , 

H ^ z ^FTim (22) 
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we rewrite Eq.(19) as 
dH 



du 



u,z) 



1 



[z + v{u)]' 



dG 
du 



(u,0) + < v(z) > 



(23) 



From the asymptotic formula Eq.(20) combined with the definition (22), it 
follows that 

lim H(u,z) = -. (24) 

Using this result we can integrate Eq.(23) over the velocity space obtaining 
an explicit expression for the mean velocity 



< viz) > 



+°° du dG ( 
-co [z + v{u)] 2 du 



u,0) 



du 



[z + is(u)f 



(25) 



The initial condition appears in Eq. (25) through the derivative (see Eq. (13)) 
dG 



du 



(u,Q) 



dcf(c, 0) 



dcf(c,0) 



(26) 



It is worth noting that we could calculate the first moment of the probability 
distribution < v(z) > without yet knowing the complete solution of the Boltz- 
mann equation. If f(u, 0) is an even function of velocity, then the derivative 
(Eq. 26) is an odd function, and from Eq. (25) we get < v(z) >= 0, and hence 
the mean velocity < v (t) >= 0, for any t > 0. 

The differential equation, Eq.(23) can be easily integrated with the use of 
the boundary condition, Eq.(24). The relation Eq. (22) between G(u,z) and 
H (u, z) yields then the complete solution of Eq.(19) 



G(u,z) = [z + u{u)\ 



dw- — — 

oo [Z + flw)\ z 



'dG 
dw 



(w,0)+ < v(z) > + 



(27) 



Differentiating Eq.(27) twice with respect to u we find the Laplace transform 
f(u, z) of the solution of the Boltzmann equation 



f(u,z) = (f)(u 
+ 



1 



dG 



z ~^ I-oo [z + v(w)} 2 \^dw 



(w, 0)+ < v(z) > 



z + v(u)} 



f(u,0). 



(2? 



Eq.(28) when integrated over the velocity yields / dvf(u,z) = 1/z, a result 
consistent with the conservation of the probability. 
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The use of the identity 

1 



+ 



dw- 



dis(w) 



[z + v(u)} J-oo [z + u(w)] 2 dw 

in Eq.(28) allows us to perform the inverse Laplace transform term by term. 
In this way we determine the structure of the time dependent solution of the 
Boltzmann equation 



/(«,i)=^) + [/(n,0)-^)]e-^ 



dw 



dG . du{w) 
dw ' dw 



te~ tu{w) + A(w,t) 



(29) 



where A(w,t) is the inverse Laplace transform of 



A(w,z)= V 30 
[z + z/(w)J 2 



Note that when the initial distribution f(u, 0) is an even function of u (and 
thus < v(z) >= 0), A(w,t) vanishes. 

The first term on the right-hand side of Eq.(29) represents the stationary solu- 
tion, whereas the second term describes the exponential decay (with relaxation 
time l/v(u)) of the initial deviation from the stationary distribution. For high 
velocities the collision frequency u(u) goes as \u\ so that the associated relax- 
ation time approaches as H" 1 . Consequently, the high energy tails of the 
distribution function converge faster towards the stationary state than the 
rest of the distribution. This phenomenon is illustrated in Fig. 4 which shows 
the evolution of the velocity distribution starting from a "gate" function. The 
initial discontinuities in velocity space are present at any finite time. Their 
position remains fixed, but their magnitude decreases with time. Finally, we 
note that the last term on the right-hand side of Eq.(29) vanishes both for 
t — > 0, and for t — > oo. 



5 Velocity autocorrelation function and diffusion coefficient 

The velocity autocorrelation function is defined by 

C[t) =< v(t)v(0) > (31) 
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Fig. 1. Velocity distribution function f(u,t) versus u at different times 
t = 0.1,l,2,5,oo. 

where the brackets denote an average over the stationary Maxwell distribution 
(Eq. 11). We can calculate C(t) exactly using the solution (28). Specifically, 

/+oo 
vF(v,t)dv (32) 
-oo 

where F(v, t) is the solution of the Boltzmann equation corresponding to the 
initial condition 

F(v, 0) = v<j>(v) (33) 

The Laplace transform of C (t) 

r+oo 

C(z) = / vF(v, z)dv =< V{z) > (34) 



is thus equal to the Laplace transform of the mean velocity < V(t) > corre- 
sponding to the solution F(v,z) (see Eq.(25)). Note that the initial condition 
(33) does not represent a probability distribution. One finds in this specific 
case the relation 

dG , , dv(u) , . , . 

ftr(«-°)=»-sr-^- (35) 
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a 



0.8 0.9 1 



Fig. 2. Reduced diffusion coefficient as a function of the coefficient of restitution: 
the horizontal line corresponds to the exact solution of the Boltzmann equation and 
the curve to the first Sonine approximation. 



Inserting Eq.(35) into Eq.(25) after integration by parts we get 

\ -i 



C(z) =< V(z) >= 2 



+oo 



du 



[z + u(u)\' 



— z. 



(36) 



Eq.(36) is an exact formula for the velocity autocorrelation function. The 
above derivation is by far simpler than that given in 0]. 

The dimensionless diffusion coefficient D* = C(0) can be now deduced by 
taking the limit z — > in Eq.(36). We find 



30 du 
o u(uY 



(37) 



Restoring the physical dimensions of the variables yields for the diffusion co- 
efficient D the formula 



D 



1 aT 



D" 



(3* 



pV M 

The numerical value of D* is 0.464139.... It is worth noting that D is propor- 
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tional to the square root of the coefficient of restitution and hence vanishes 
when a — * 0. 

The diffusion coefficient for a tracer particle has also been obtained approxi- 
mately from a normal Chapmann-Enskog solution of the Boltzmann equation 
H, k| yielding 



D* = v — -. (39) 

3 + 2a - a 2 K ' 

Figure 2 shows the reduced (dimensionless) diffusion coefficient D* versus the 
coefficient of restitution a. 

We note that the diffusion coefficient for elastic (a = 1) hard rods in one- 



dimension is known exactly: D* = \/\/2n ~ 0.3989. .|14j. Probably fortu- 
itously, at a = 1 the lowest order Sonine approximation to the solution of 
the Boltzmann equation gives a result for D* that is closer to this exact re- 
sult than the prediction of the Boltzmann equation itself. When a < 1, the 
reduced diffusion coefficient within the Sonine approximation increases when 
a decreases, whereas the result (Eq. 37) provided by the Boltzmann equation 
is independent of a. The two curves cross at a = 0.574. In both cases, of 
course, the full diffusion coefficient D vanishes when a goes to zero due to the 
square- root dependence on the coefficient of restitution in Eq.(38). 



6 Simulation results 



In order to study the situations where M/m ^ ct, we have solved the Boltz- 
mann equation numerically by using a DSMC method 0|. We focus here on 
the stationary dynamics. Figure 3 displays both the numerical results of the 
simulation and the exact results obtained by an inverse Laplace transform of 
Eq.(36): the ratio of the tracer particle over the mass of the bath particle is 
equal to a, when a = 0.6. Note that decreasing the coefficient of restitution 
increases the relaxation time of the normalized velocity autocorrelation func- 
tion but does not significantly change its shape. The dashed curve corresponds 
to the exact solution, Eq. (36), and matches very accurately the simulation 
results when a = 0.6. 

As a decreases, the velocity autocorrelation function decreases more slowly. 
This behavior is, at first glance, counterintuitive. When a is small the postcol- 
lisional velocity of the tracer is smaller than when a is larger. One therefore 
expects that the correlation function should approach zero more rapidly when 
a is small. This reasoning, however, does not account for the time between 
collisions which increases as a decreases. Evidently it is this effect that domi- 
nates, leading to the observed behavior. 
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Fig. 3. Normalized autocorrelation function C{t) of the velocity versus time 
for M = 0.6m and for various values of the coefficient of restitution 
a = 1.0,0.8,0.6,0.4,0.2,0. The dashed curve corresponds to the inverse Laplace 
transform of Eq.(36) 

7 Self correlation function 

The more general case of spatially inhomogeneous states also admits an ana- 
lytic solution of the Boltzmann equation. Since the strategy is basically similar 
to that developed in Section 3, details of calculations are given in Appendix 
A. The final result for f(q,u,z) reads 





1 



1 



+ 





1 (dG_ 
[z + u(w) + iqw] 2 \dw 



{q,w,Q) 



dw 



(w)+ < v(q, z) > 



)] 

(40) 



Here G denotes the Fourier transform 




(41) 



12 



and G(r,u,t) is denned by Eq.(A.l) in analogy with the definition (13). 



The conditional density F(r, t) can be obtained from the inhomogeneous Boltz- 
mann equation by integrating f(r,v,t) over the velocities assuming that 



f(r,v,0) = 5(r)<j> M (u) 



(42) 



i.e. that initially the tracer particle was located at the origin with the equilib- 
rium velocity distribution. This is in fact the definition of the self-correlation 
function 

/+oo 
dvf(r,v,t). (43) 
-oo 

As, in accordance with the definition (A.l), 



2f(r,v,t) = ^G(r,v,t) 



(44) 



one finds that the Fourier-Laplace transform of F(r,t) is given by 
F(q,z) = j dre-^ jf dt e~ zt F(r,t) = - —(q j+ oo,z)-—(q,-oo,z) 



From Eq.(A.13) a simple expression follows 

~ 1 -iq < v(q,z) > 
F{q, z) = , 

whereas using the initial condition Eq.(42) one finds the formula 

~ z(l+q 2 )J(q,z)-2 



< v(q, z) >= iq 



with 



Finally, one gets 



J(q, z) 



z(l + q 2 )J(q, z) + 2q 2 

fo ° dw 

oo [z + iqw + u{w)] 2 

2(1 + q 2 ) 



F(q,z) = i + £ (1 

V ^ z z \ z(l + q 2 )J(q,z) + 2q 2 / 



(45) 
(46) 

(47) 
(48) 

(49) 



The wave vector dependent poles of F(q, z) satisfy the equation 

z(l + q 2 )J(q,z) + 2q 2 = (50) 

For the hydrodynamic pole z(q), which tends to zero when q — > 0, Eq.(50) to 
dominant order yields 



zJ(0,0) + 2q 2 = 
z + Dq 2 = 



(51) 
(52) 
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where D is the diffusion coefficient derived in section 5 with the use of the 
Green-Kubo autocorrelation formula (see Eq.(37) ). Here it appears in the 
hydrodynamic long-time and long-length scale limit of the solution of the 
Boltzmann equation which expresses the density fluctuation of the tracer par- 
ticle. 



8 Conclusion 



We have shown that, within the framework of the one-dimensional Boltzmann 
equation, the dynamics of a granular tracer particle can be mapped onto the 
dynamics of an elastic tracer particle with a renormalized mass. A recent 
article of Santos and Dufty [16[ shows that this result can be extended to 
more general situations. 

J. P. acknowledges financial support by the CNRS, France, and the hospitality 
at the Laboratoire de Physique Theorique de la Matiere Condensee, UPMC 
(Paris) where this research has been carried out. 



A Inhomogeneous solutions of the Boltzmann equation 



We summarize the basic steps leading to the solution of the inhomogeneous 
Boltzmann equation. The quantity which generalizes Eq.(13) is the function 
G(r, u, t) defined as 



G(r, u, t) 



+00 



dc\u — c\f(r, c, t). 



(A.1) 



Let us denote G(q,u,z) the Fourier-Laplace transform of G(r,u,t), 



G{q,u,z) 



+00 



dre iqr 



+00 



dte- zt G(r,u,t). 



(A.2) 



The Boltzmann equation, Eq.(5), can be then rewritten as 

r • / 2 G . . d 2 v(u)~, . d 2 G 

[z + iqu + ,/(«)] — (g, u , z) - -^rG{q, u, z) = — (g, u, 0) 

Integrating Eq.(A.3) over the velocity interval ] — 00, u] yields 



(A.3) 



[z + iqu + z/(u)]^(g, u, z) - 



dv{u) 
du 



+ iq 



G(q,u,z) 



dG 
du 



(q,u,0) + B(q, z) 
(A.4) 
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where the function B(q,z) is to be determined. Similarly to Eq.(22), it is 
convenient to introduce the function H(q,u, z) 



H(q,u,z) 



G{q,u,z) 



[z + iqu + v(u)] ' 



(A.5) 



which satisfies a differential equation analogous to Eq.(23) with u(u) replaced 
by \y(u) + iqu], and < v(z) > replaced by < v(q, z) >. Therefore, the solution 
of Eq.(A.3) can be expressed as 



G(q, u, z) — [z + v{u) + iqu] 



(A.6) 



dw 



'dG 



oo [z + v{w) + iqw] 2 \dw 



(q,w,0) + B(q,z))+D(q,z) 



The functions B(q,z) and D(q,z) can be determined from the boundary 
conditions satisfied by G(q,u,z) and H(q,u,z). The asymptotic behavior of 
G(q, u, z) for large u at fixed z and q is given by 



G(q, u, z) ~ \u\ 



+00 



dwf(q, w, z) — sgn(u) < v(q, z) > 



where 



< v(q, z) > 



+00 



dwwf(q, w, z) 



(A.7) 



(A.8) 



The continuity equation (conservation of the probability) is expressed as 

r+00 _ 

dwf(q, w, z) +iq < v(q, z) >= 1 (A. 9) 

By combining Eqs.(A.7) and (A. 9), one obtains the asymptotic relation 



G(q, u,z) ~ — [1 — iq < v(q, z) >] — sgn(u) < v(q, z) > (A. 10) 



The limit of H (q, u, z) when u — > ±00 is then given by 

l — iq<v(q,z)> 
hm H(q,u,z) = J }\' ' A.ll 

By integrating the differential equation obeyed by H (q, u, z) over velocity, one 
obtains a sum rule analogous to Eq.(25): 

-/ \ / 2ig r+°° du dG , n \ . . 

«*'^>=-{TT? + L lz + V (u) + , q u]^ ' U ' 0) ) (A - 12) 



du 2q 
+ 



2 \-i 



[z + u{u) + iqu] 2 1 + q- 
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The knowledge of the boundary conditions (A. 10) and (A. 11) permits one to 
write the solution (A. 6) in an explicit form 

G(q, u, z) = [z + u{u) + iqu] (A. 13) 

r u , 1 (dG, . . \ 1 - iq < via, z) > 

x / dw \ i i rr — 12 o-(?. u; > )+< u > + — n — -\ 

7-oo [z + u(w) + iqwf \ow J {l — zq)z 

Taking the second derivative of Eq.(A.13) with respect to the velocity variable 
u leads to Eq.(40). 
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